A Superconducting instability in the surface of a topological insulator 
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It is shown that a superconducting instability appears in the electronic states on the surface of a 
topological insulator due purely to electromagnetic interactions. The discussion of this instability 
is based on the analysis of the gap equation using the repulsive term coming from the axion term 
and the assumption that the system behaves like a two dimensional Fermi liquid. It is shown 
that this superconducting instability appears in the p-wave channel and an estimate of the critical 
temperature of the transition is given. 



Introduction. The topological insulating (TI) state of- 
fers the possibility to study physical phenomena that 
were considered as exclusive of the High Energy Physics 
not so long ago[ll-0]. The prominent example is the 
topological magnetoelectric coupling known in the former 
context as the axion term 0, . This term has proved 
to be quite attractive to the Condensed Matter commu- 
nity due to its fascinating properties like the possibility of 
the existence of Majorana Fermions in superconductor- 
TI junctions and a repulsive Casimir interaction between 
topological insulating plates d, Q (or an up-to-date re- 
view of the subject, see Q). Another interesting exam- 
ple, a superconducting state on the surface of a three 
dimensional TI, has been analyzed in a context of a gen- 
eral symmetry analysis or in basis of a generic short range 
interactions (9I [loj . From the experimental side, a super- 
conducting phase has been observed Bi2Se3 doped with 
Cu[lllll2j. It is observed that the presence of a Cu layer 
is the key ingredient to induce superconductivity. Here 
our aim is not to explain these experimental observations 
but rather suggest a microscopical mechanism of super- 
conductivity when the TI surface states are magnetically 
gaped being the magnetoelectric term the ultimate re- 
sponsible for the electronic pairing. The problem is of 
considerable interest for several reasons: it constitutes 
a physical realization of the residual pairing mechanism 
proposed to occur in the context of the Quantum Hall 
Effect (QHE) at half filling [l3l|. and also the resulting su- 
perconducting order parameter enjoys p-wave symmetry 
opening the possibility of finding majorana states at the 
vortices in this system, in this case without relying in any 
proximity effect jg, [3] . 

The model. Our starting point is the action for the 
electronic states coupled to the electromagnetic field: 



S = J {S{x3)Csurf + Cem 



-'bulk 



(1) 



We will assume that while the surface electronic states 
live on the surface of the TI (this is the origin of the term 
S{x3)) both the bulk gapped electronic states and the 
electromagnetic field live in the three spatial dimensions. 
This fact implies that we can not longer describe the elec- 
trodynamics of the surface states by a proper QED3, as 



is done in |15l . Ila | and some crucial differences arise when 
the photon propagator is calculated, as we will see. An- 
other important difference is that although the spectrum 
of the electronic states is well described by a Dirac equa- 
tion, the typical Fermi velocity of these states is around 
c/600[17j and then we expect that retardation effects are 
negligible. We will use (3 = vp/c as a perturbation pa- 
rameter to keep terms that will lead to the effective in- 
stantaneous Coulomb potential. 

First of all, we will focus on the electromagnetic part 
of the action in ([1]). Because we are mainly interested 
in the non retarded effects of the Coulomb interaction in 
the surface states, we will not explicitly write the term 
corresponding to the gapped bulk states. Instead of that, 
we will integrate out the bulk fermions taking them into 
account in the modification induced in the action for the 
electromagnetic action. It is already well known that 
this contribution is twofold: on the one hand, we con- 
sider the system to be a nonmagnetic insulator (/i ~ 1), 
that is, the effects of the bulk electronic structure (and 
other effects at higher frequency scales) are encoded in 
the dielectric susceptibility e(a;,k). Moreover, here we 
are dealing with physical effects at a very low energy 
scale, smaller than the gap energy of the bulk electrons 
and other polarization effects so it is reasonable to ap- 
proximate e by its value at zero frequencies and momenta: 
e(a;,k) ^ e(0,0) = e^. On the other hand, an axion- 
like term also arises in the electromagnetic response of 
the bulk electrons in a Tl[4|. In order to make the ax- 
ion term physically observable we need to introduce a 
time reversal breaking element in the surface of the TI, 
like a magnetic field, a ferromagnetic coating or mag- 
netic impurities. The side effect of such an element is to 
open a gap in the surface electronic spectrumjjj]. This 
term, called magnetoelectric term, strongly modifies the 
Maxwell equations in the material medium [5'] and thus 
it becomes evident that if we want to consider the elec- 
tromagnetic field propagating in the three spatial dimen- 
sions we must take into account the presence of such a 
term and thus it will modify the long range Coulomb 
interaction among electrons. 

We will follow the strategy settled in [3] to compute 
the EM propagator. Because the electromagnetic field is 
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coupled to a a current constrained to live on a surface, 
we can get rid of the fourth component of the electro- 
magnetic field Afj^{x) by using a suitable choice of gauge. 
In this context we choose the Feynman gauge d^A^j^ = 
because this particular choice of gauge not only allows 
us to decouple the dynamical evolution of A3 from the 
other components but it makes the even part of the pho- 
ton propagator to be proportional to the surface metric 
tensor rjij = diag(\, —1, —1). 

Before discussing the action of the electromagnetic 
field Aj{x) we can take some advantage of the fact that 
we are considering a constant dielectric function tr- We 
anticipate that the effective Coulomb interaction at low- 
est order in perturbation theory in the charge e will be 
of the order of . Also we know that the dielectric con- 
stant always appears in the form of the ratio /tr = otr 
so we can safely work with a free electromagnetic action 
throughout the intermediate calculations and then just 
substitute in the final result the bare charge by this ra- 
tio. 

Under these considerations, the effective action for the 
electromagnetic field including the axion term will be: 



Sem = j d^x {A^ri'^'d^A, + aeix)e''''P^d^A,dpAx) , 

(2) 

where a is the fine structure constant and 0{x) is the 
axion field. We can model the axion field for a perfect 
interface vacuum- TI with a step function in the z di- 
rection: d{x) — 6Q{x3) where now 9 takes the value of 
(2n -I- 1)tt, n £ N. Integrating by parts the second term 
in ([2|) and integrating out A3, Sem becomes 



S 



EM 



d'^xiArj'^d^A, - 5(x3)n3e'*-''A,9,Afc). 



2tt 



(3) 

In ([3]) we see that the axion term acquires the form of a 
Chern-Simons (CS) termj^j. However the dependence of 
the electromagnetic field of the coordinate 2:3 radically al- 
ters the form of the electromagnetic propagator Dij {x,x') 
and it will be completely different to what one expects 
from the electromagnetic propagator in a CS QED3[l^. 
The first thing we have to notice is that now the system 
is no longer translationally invariant along the direction 
X3 so the spatial dependence of the propagator will be 
Dij (x — x' , a;3 , 0:3 ) . Because now the electromagnetic field 
is coupled to a system which lies on a surface, we do not 
need to know the complete dependence of Dij with X3 
and X3 but only the value at the surface of Dij{k, 0,0) 
(the Fourier transform in the rest of the coordinates). 
By noticing that the second term in is proportional 
to S{x3) and using the Dyson equation for the electro- 
magnetic propagator we get: 



Aj(g,0,0) 



^0 



^\q,0)+i^e^i^qi' 



(4) 



In the last expression, D^^~^{q,0) = Dl^~^{q,X3 — x'^ = 
0) stands for the integral of the free photon propagator 
which in momentum representation reads(19i]: 



(5) 



and \q\ = ^/qq — q^. Note that now both terms in the 
denominator of (jlj are linear in q. Then using (|4]) and 
([5]) we arrive at the expression for Dij {q, 0,0) that we 
will use to compute the effective Coulomb interaction (in 
the nonretarded limit qo 0): 



A,(q,0,0) 



2(1 + 



Vij 



iaO q 



(6) 



The effective propagator for the electromagnetic interac- 
tion is not equal to the propagator for a purely CS QED3 
theory. In the present case, the effective propagator is 
still massless, despite of the presence of a CS term, con- 
trary to what happens when the photon only propagates 
in two spatial dimensions, where the CS term induces a 
mass for the photon. 

Let us turn our attention to the action of the sur- 
face electronic states. In order to make the axionic term 
observable we have to break the time reversal symme- 
try in the system. This situation can be experimentally 
achieved by directly adding a magnetic coating on the 
surface or by adding enough magnetic impurities in the 
bulk of the Tl[3|. The side effect is that we induce a 
mass in the surface electronic spectrum. In the action 
([T]) we will define a;° — ct so that the action term for the 
surface states is 

Ssurf = / d'*a;(5(x3)?/5 {-i-/iAP^{dj + ieAj) - m) i/). 



(7) 

In ([7]) the indices {i,j) run only from to 2 and we have 
defined the matrix M*^ = diag{l, —f3, —/3). Since the 
parameter /3 is already defined, we will set c = 1 in the 
rest of the paper. The matrix convention used in ([7|) is 
7o = o'2,7i = iffy, and 72 = —iux, so the spin operator 
for the surface estates is S = z x a ^ —ij. Going to a 
second quantization formalism, we can expand the field 
operators in terms of creation and annihilation operators 
of states in the valence (6k) and conduction bands (ok): 



V^(x) 
^(x) 



E 

k 

E 



-zkx 
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(8) 
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where Uk and are the planewave solutions of the Dirac 
equation resulting from (O (wk = 7iMk)- 



part of the effective propagator, the electron wavefunc- 
tions Uk and the matrix M*-' : 



Wk = 



j m + uji, / /3(fcx - iky) ^\ 



In what follows we will set the Fermi level slightly above 
the bottom of the conduction band so in virtue of the 
smallness of /3 we can write Wk — /i = -yZ/^-^k^ + m'^ — m w 

/32k2/2TO = Ek. 

According to the exchange of a photon induces 
an effective electron-electron interaction containing two 
terms corresponding to the even (proportional to rjij ) and 
odd (proportional to eiji) terms in the effective propa- 
gator. The even part gives rise to the standard non- 
retarded Coulomb interaction. As it happens in stan- 
dard two dimensional electron gases the existence of a 
finite Fermi surface will induce a polarization screening 
of the Coulomb interaction and hence a quasiparticle life- 
time proportional to lo'^Ihuj, for energies u close enough 
to the Fermi level, so we can safely consider a Fermi liq- 
uid picture where the Coulomb repulsion among electrons 
leads to a weakly interacting system of fermionic quasi- 
particles with renormalized mass m and Fermi velocity 
/3 and where an extra interaction term, corresponding to 
the odd part of the effective propagator ([5]) remains. Of 
course, if we insist in the static screening of the Coulomb 
interaction we are forced to accept that the odd part of 
the propagator also changes by polarization effects so we 
need to know how the bare odd term gets modified by 
the inclusion of the polarization term ny(a;,q) in ^ 
in the static, long wavelength limit. Both the even and 
odd parts of the polarization function Ilij{uj,q) are well 
known in the literature j20|. Gauge invariance tell us 
that the polarization tensor is defined through two scalar 
functions, n'^"^"(w,q) and 11°'^'^ {w , q) . At zero frequency 
and in the long wavelength limit they take the values 
n™™(0,q ^ 0) ~ 2^i/^^ = A (A is the Thomas-Fermi 
wavevector) and n°''''(0,q ^> 0) ~ m/^ so we can write 
the following expression for the inverse of the odd part 
of the effective propagator (7 = ar{0 + er^)/27r): 



(10) 



Because now we are interested in the regime of energies 
close to the Fermi level, we will only take into account 
terms in the interaction hamiltonian corresponding to in- 
traband processes, so the relevant interaction term com- 
ing from the odd part of the dressed effective propagator 
reads in the Cooper channel q = k' — k: 



Hint = y^F(k,k')a+ak'alkQ-k', (H) 

k,k' 

where the potential ^ (k, k') is defined through the odd 



T/(k,k') ^ arD°f'^{0,k' - k)M*''A/p%^7^Ukmlk7.W-^k'- 

(12) 

In order to find the form of l^(k, k') we need to know 
the form of the vector ri(k, k') = wJ^7iWk'. Using the 
explicit form of Uk in (O up to order we have 

ro(k, k') = 1 - 1^ (|k - kf + zz(k' X k)) , 
ri(k, k') = ^{k' + k*), r2(k, k') = ^(fc' - r ),(i3) 

2m 2m 

where k = kx + iky and k* = kx — iky. It is worth to 
mention that the zero component of F.^ is even under the 
simultaneous change k, k' — k, — k' while the spatial 
components 2) are odd. This fact prevents the sum- 
mation in (fT2|) to be zero. 

With all these considerations in mind, we arrive to the 
explicit form of the effective potential ([l] 



ar7/3^ 



|k-k' 



2iz(k' X k) 



m (|k-k'| A)2 -|-72|k-k'|' 



(14) 



The equation for the gap Ak can be easily written by 
using an standard Bogoliubov transformation and taking 
into account (fTT|) and : 



<Pk' 



A. 



47r2 



= F(k,k') 



(15) 



Results. The presence of the complex unity in the sec- 
ond term in ()14p prevents Ak to enjoy s-wave symme- 
try. This is expected, however because the time reversal 
symmetry is broken due to the presence of the magnetic 
coating, responsible of the mass of the surface electrons. 
We will employ a p-wave symmetric ansatz for the gap 
parameter: Ak — A{kx + iky) = A|k|e''^''. Also, the 
Thomas-Fermi wavevector A introduces a characteristic 
lengthscale in the problem. Because of that, two differ- 
ent regimes in eq. (jl5p allows for analytical tractability: 
k|, |k'| ^ A and |k|, |k'| » A, leading to upper and lower 
bounds of the transition temperature Tc respectively. 

In the regime |k|, |k'| ^ A we can substitute the de- 
nominator in (|14p simply by A^. By the shift 0k' 
'/'k' — 0k the angular integration in (|15p is straightfor- 
ward leading to: 



A 



8vr2A2 Jo 



dlk'l 



Alk' 



(16) 



Because we are formally working in the limit A — 00 we 
can substitute the upper limit of the |k| integral A by the 
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Thomas-Fermi wavevector A and keep the leading term 
in the limit /3 ^ 0. The result is simply 



A = 



67r 



(17) 



In the most favorable situation where (3 = 1/600, /i ~ 
m « 7meV, 9 = n, and e,, 1, the critical temperature 
is estimated in this limit to be Tc ~ -A « 6K. 

KB 

In the opposite regime |k|, |k'| » A, the lengthscale A 
disappears from (jl4p . As we have mentioned in the intro- 
duction it is interesting to note that such limit strongly 
resembles to the attractive residual interaction among 
quasiparticles in the half filling regime in the QHEfls"]. 
After performing the angular integration, eq. (jl5p takes 
the form 



Alkl = 



|k| 



d|k'| 



Alk'l 



Ik'l 



A2 + /^7|k'|2 

4m^ I I 



+ 



+ / d\k'\ 
/|k| 



Alk'l 



A2 + /i^|k'|2 

4m^ I I 



(18) 



The second term in the right hand side is logarithmically 
divergent and dominates the whole right part of (1181) . 
Using a hard cutoff A, the leading contribution to A is 
approximately independent of |k|: 



(19) 



In the last expression the parameter 7 is controlled by 
the value of a^, so when 1, 7 ~ a ~ 1/137 (now the 
most unfavorable situation) the exponential is extremely 
small, meaning that the lower bound for Tc is nearly zero. 
However, because we are working with a long wavelength 
effective theory the dominant contribution to the super- 
conducting order parameter Ak comes from small ener- 
gies around the Fermi level, implying that a reliable es- 
timation for Tc is the one found in the small momentum 
regime. We have to mention that the p-wave symmetry 
of the order parameter is consistent with two basic facts: 
as we said, the time reversal symmetry is already bro- 
ken because the magnetic coating and because the spin 
polarization of all electrons is the same, determined by 
the sign of m, so the origin of the p-wave symmetry of 
Ak and the antisymmetry of the Cooper pair is purely 
orbital. 

Because of the p-wave symmetry of the superconduct- 
ing order parameter and the dimensionality of the prob- 
lem, we can speculate with the fact that the supercon- 
ducting transition will be of Kosterlitz-Thouless (KT) 
type, in a similar way to the transition found in super- 
conducting thin films [21I, [S^]. The possible existence of 
vortices in A(r) might imply the existence of Majorana 
fermion states at the core of such vortices, according to 
ref. 23|. This opens another route for searching Majorana 
zero modes in topological insulators. 



Conclusions. It is found that the topological magneto- 
electric term present in a three dimensional topological 
insulator strongly modifies the electron-electron interac- 
tion when it becomes observable. The Induced Chern- 
Simons term gives rise to an attractive term that eventu- 
ally induces a superconducting instability in the surface 
electronic states when the standard Coulomb repulsion is 
(statically) screened. The superconducting order param- 
eter enjoys p-wave symmetry, which is consistent with the 
fact that time reversal symmetry is broken and the an- 
tisymmetry of the pair wavefunction is not determined 
by the spin part of the electronic wavefunctions but it 
has orbital origin. We have found an upper bound for 
the transition temperature Tc to be of the order of 6K. 
Finally, we have commented about the possibility of KT 
character of the transition and the existence of Majorana 
zero modes at the core of the vortices appearing in A(r). 
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